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Abstract 

This work is devoted to the study of a compressible viscoelastic fluids satisfying the 
Oldroyd-B model in a regular bounded domain. We prove the local existence of solutions 
and uniqueness of flows by a classical fixed point argument. 

1 Introduction 

In this paper, we study the local existence of solutions for compressible viscoelastic fluid flows 
in the case of the Oldroyd-B model in a regular bounded domain in M 3 . We also show the 
uniqueness of solutions. We prove l'existence by using the classical method based on the 
Schauder fixed-point theorem. Valli, in [8], show the local existence in the case of the Navier- 
Stokes equations. The case of the Oldroyd model for incompressible fluid is studied by Guillope 
and Saut in [3]. Talhouk shows the existence and the uniqueness for Jeffreys model's in [6|. 

This paper is organized as follows. Section |2] is devoted to the modeling of the problem 
and to the definition of well-prepared initial conditions. The principal notation and results are 
detailed in Section [3j The local existence of regular solutions is given in Section 4. 

2 The Modeling 

2.1 Unsteady Flows of Compressible Viscoelastic Fluids 

Consider unsteady flows of viscoelastic fluids in a bounded domain S7* of M. 3 with a regular 
boundary T* . The system, obtained from the laws of conservation of momentum, and of mass, 
and from the constitutive equation of the fluid, reads as follows jl]: in Q^* = (0,T*) x fi*, 

P* (j^r + (u* • V*)u*) = p*f*+div*(r*-p*I), 

^- + div*(p*u*) = 0, (2.1) 



T>t* V ^ Vt* J 
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The *-variables are the dimensional ones in the domain of the flow Q*, and T* > is a 
dimensional time. The unknowns are the velocities u*, the density p*, and the symmetric 
tensor of constraints r*. r] is the total viscosity of the fluid, A > is the relaxation time, and 
/i is the retardation time (0 < jx < A). 
V a r* 



is an objective derivative of the tensor r*, given by 

+ (u* • V*) ) t* + t*W* - W*t* - a(D*r* + r*D*), 



Vt 

V a r* ( 8 



Vt* V dt 

where W* = W*[u*] = ^(V*u* - V* T u*) and D* = D*[u*] = ^(V*u* + V* T u*) are, re- 
spectively, the rate of rotation and the rate of deformation tensors, a is a real parameter in 
[-1,1]. 

System (12. ip is completed by a condition on the boundary, 

u* = on S^* = (0,T*) x r*, 

and by the initial data 

u*(0,-) = u$, p*(0,-)=p* , T *(0,-) = r *, inO*. 

We split t* into two parts: the Newtonian one r* related to the solvent, and the polymeric 
one r*. We may write 

t =r s +r p = 2r] s T> + r e , 



where r* = r* — ^^-div*u*j I, and I is the identity tensor. rj s = rjp/X and £ s are the solvent 
viscosity and the group viscosity, respectively. Since we are interested in a model for weakly 
compressible fluids, we suppose that £ s = 0. From the third equation in (|2.1|) . we can deduce 
that t* satisfies the equation 

V T* 

where r/ e = ij — ij s is called the polymer viscosity. r] s and r] e are two non-negative numbers. 
Therefore, under the assumption £ s = 0, System (12. ip is equivalent to the system in Q^* , 

P* (j^T + ( u * • V*)u*^) = p*f* + Vs (A*u* + V*div*u*) - Vy + div*r*, 

^+div*(p*u*) = 0, (2-2) 
dt* yH ' 

r* + A-^— = 2ry e D*u*, 
Vt* 

where we have denoted r* by t* to simplify the notation. 
2.2 Well-Prepared Initial Conditions 

We first define the Mach number e as being the ratio of the typical velocity of the fluid U to 

/ , * \l/2 

the speed of sound ( -^{p* Q ) J in the same fluid at the same state. We divide the density 
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p* = p* e into two parts: a constant one independent of e, and a remainder, which is small 
for small e's, say 

P = A) + °( e ) = Po + e o- . 

We also suppose that the initial conditions p^ £ , u* £ and t^ £ are well-prepared, which means 
that they take a similar form, say 



a* 6 



U 



7 f 0+ O(e 2 ) = pl + s 2 <, 
v* + v* e , with divv* = 0, 

S* I Q*£ 
' 



where v* and S*, are, respectively, a vector and a symmetric tensor, both independent of e. 
We assume 

m* = min p* > and 9Jt* = max p* . 

n* n* 

Assuming that p* = p*(p*) is regular, say class C 3 at least, we remark 

/_* . 2 * \ dp* . 2 f 1 d 2 p* /_* . 



1 se a ) as. 



dp* ^* ^ ^ 

We introduce the function w*, defined by w*(a*) = ^*(p*, + e 2 a*) - ^(p*,), and remark that 
w* depends on e, satisfies w*(0) = 0, and is of class C 2 at least. 
Replacing p* by its value in the first equation (|2.2p . one infers 



-2„*\ dp* 



du 



(P* + e 2 a*) [ — + (u* • V>* ) + e'^M + e 2 a*)V*a 



■dp* 



dp* 



+eV)f* +t/ s (A*u* + V*div*u*) + div*r*. 



We can also rewrite this equality, by taking into account the definitions of w*(a*) and of the 
Mach number e, in the form 



[Po+e a 



du* 
~dt* 



+ (u* • V*)u* + (U ) 2 VV 



= (p* + e 2 a*)i* + 7] s {A*u* + V*div*u*) + dW*r* - e 2 w*{a*)V*a* 
From the second equation in (12. 2h we easily deduce 

,da* 



dt* 



+ p*div*u* + e 2 div*(<7*u*) = 0. 



Finally, System (|2,2p can be written as follows, in Q^*, 

' (p* + eV*) + (u* • V*)u*) + (U ) 2 VV* = (p* Q +e 2 a*)f* + div*r* 



dt* 



,da* 



+n s (A*u* + V*div*u*) - e 2 w* (a*)V*a* , 



e 2 — + p*div*u* +£ 2 div*(a*u*) = 0, 
T* + A-^- = 2f7eD*[u*l. 



(2.3) 
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2.3 Dimensionless Variables 



We introduce the dimensionless variables, 



x* = L x, t* = ^j-t, p* = a p, w*(a*) = (\J ) 2 w(a), 



u 



U u, a* = a a, t* = T r, p* (p*) = T p(p), f* 



(Uo) S 

^0 



and T = ^ 



where L represents a typical length of the flow. The real numbers a 
characterize the density and the stress tensor of the fluid. 0, denotes the non-dimensional 
domain of the flow, with boundary T, and T > a non-dimensional time. 

We introduce three non-dimensional numbers: a number a similar to the Reynolds number 
for incompressible flows, the Weissenberg number We, and a number u relative to the viscosities 
of the fluid, 

a = tL = mh, We = ^, 

a rj L 

We also define 

dp , 2 \ dp 



oj = 1 . 

1} 



w(a) = a < — (a + e 2 a) 
dp 



dp 



[a 



In dimensionless variables, System (|2.3p takes the form, in Qt = (0, T) x Q, 



u' + (u • V)u + 



1 



a + e 2 a 



-Va 



a' + e 2 a div u + div (cru) 



1 — w , a „ \ div r 
f H : — — ( Au + Vdiv u) + 



a + e 2 a 



e 2 w(a)Va 
a + e 2 a 



a + e 2 a 



0. 



, r + We{r' + (u • V)r + g(Vu,r)} = 2wD[u], 



(2.4) 



•ii ■ . du . da dr 

with the notation u = — — , a = —— and r = — — , and 

dt dt &t' 



g (Vu,r) =rW[u] - W[u]r - a(^D[u]r + rD[u]J . 

Introducing the differential operator A = — (A+ Vdiv) we may rewrite System (|2.4p as follows, 
in Q T , 

a(u' + (u • V)u^+(1 - u>)Au + Vex = F(u, a, r) + divr, 



a' + (u • V)cr + cdiv u = — e 2 adivu, 
r + We(V + (u- V)r + g(Vu,r)) = 2wD[u], 



with 



F(u,ct,t) 



_ (1-W)£V, 

af + Au + 

a + e 2 cr 



e 2 (a 



w(a)) 



a + e 2 o" 



e 2 a 



a + e 2 cr 



-div r. 



(2.5) 



(2.6) 
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System (|2,5p is completed by an homogeneous condition on the boundary, 



u = on E T = (0, T) x T, 



(2.7) 



and by three initial conditions 



u(0,-) 



U , cr(0,-)=0"o, t(0, -)=t , 



in O. 



(2.8) 



We also assume the followings 



0<m 1 = — <a + eV < SDTi = 



in fi, 




where tri! and dJt 1 are some given constants. 

3 The Notation and Main Results 
3.1 Notation 

f2 is a bounded domain in IR 3 , with a regular boundary T, and n denotes the unit outward- 
pointing normal vector to T. For x = (x 1 ,^ 2 ,^ 3 ) G M 3 , we denote by |x| its Euclidean norm. 

We will use the following spaces: the Lebesgue spaces L p (f2), 1 < p < +oo, with norms ||-|| LP 
(except for the L 2 (f2)-norm, which is denoted by ||-||); the Sobolev space H fc (f2), k G N*, with 
norm ||-|| fc and inner product ((•, •))&; the vector spaces L 2 (f2) and H fc (0) of vector-valued or 
tensor- valued functions with components in L 2 (f2) and H fc (f2) respectively, their norms being 
denoted in the same way as above. We will also use the homogeneous Sobolev space Hq(O) and 
its dual H -1 («). 

If / is an interval of R+ and k G N, C(/;H fc (ft)) is the space of vector- or tensor-valued 
functions which are continuous on / with values in H k (Q). The norm, in this space, is denoted 
by Ct(/;H fc (r2)) is the space of functions of C(/;H fc (f2)) which are bounded on /. 

The space L p (/; H k (ft)), for 1 < p < +oo, and k £ N, consists of p-integrable functions on / 
with values in H fc (^)- For 1 < p < +oo, k G N and < T < oo, the norm in LP((0, T), H fe (fi)) 
is denoted by [ • ] P) k,T- Lfoc^+J H fc (fi)) is the set of functions which are in L 2 (I; H. k (fl)) for all 
bounded interval / in M + . 

The letters C, c% or q, i, j = 1, 2, • • • , will denote constants taking different values, but not 
depending on e. C n will be a constant, taking different values, and depending only on f2. (2.1) n 
denotes the n— th equation of System (2.1). 



3.2 The Main Result 



Recall the problem under study: 



< 




F(u,<t,t) + divr, 
—e 2 a div u, 
2wD[u], in Qt, 



(3.1) 



u(0, •) = u , 
cr(0, •) = cr , 
r(0,-) = r , 



in 0, 
in f2, 
in f2, 
on Ex 



u = 
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where F is defined by (|2,6p . 

Theorem 3.1. (Existence of a local solution) Assume £1 C M 3 is a domain of class C 3 . Let m 1 
and 9Jti be two real constants such that < m 1 < ^Xft 1 . Assume 

f G Lj^M+jH 1 ^)), with? G Lf oc (IR + ;H- 1 (fi)), u G H 2 (Q) n Hjjj(fi), r G H 2 (0), 

(To G H 2 (f2), u»ih / cro(x)dx = 0, and < m 1 < a + E 2 ao < DJli , in f2. 
Jn 

Then there exists a time T 1 > and a solution (u, a, r) of Problem K2. 5] )- ( TO)) in Qti = 
(0,T\) x f2, satisfying 

u G L 2 (0, T i; H 3 (0)) n C([0, Ti]; H 2 (ft) n HJ(fi)), 
u' G L 2 (0,T i; HJ(O)) flCdOJ^L^fi)), 

(t,<t) G CdCT^H 2 ^) x H 2 (ft)), (r>') G C([0, T,]; EP(fi) x H 1 ^)), 

with 

/ <r(-,x)dx = 0, m [0, Ti], and <a +e 2 a <2dJl 1 , in Q T . 
Jn 2 

Theorem 3.2. (Uniqueness of a local solution) There exist a unique solution of Problem \2. 5\) - 
\2.8\) . given in Theorem \3.1\ 

To show that the local solution found in Theorem 13.11 exists for all times under certain 
regularity and smallness conditions on the data, we also assume that the function w G C 2 (R) 
has the following properties: for all h G L 2 (0, T; H 2 ($7)), 

< C||fc'|| , | KM 1 1 < C INI > 

\\ W {h)\\ k <C\\h\\ k , k = l,2, (3.2) 

for some constant C depending on £1 and w. 

Remark 3.3. There are several examples of functions p = p(p), for which w satisfies the 
conditions above. Let us quote the case where the pressure is given by the linear state law 
p(p) = \{p — a), as well as the case of isothermal compressible perfect fluids, where p(p) = 
(C s ) 2 p, and C s is the velocity of sound in the fluid. 

4 Existence and Uniqueness of Local Solutions 

We prove Theorem 13.11 by using the classical method based on the Schauder fixed-point theo- 
rem. To do that in our case, we study three linear problems: the first one has the velocity u 
as unknown, and the next ones are two transport equations for the density a and for the stress 
tensor r respectively. The parameter e is fixed in the interval (0,1]. 
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Let w, 7r and ip a given vector, function and the symmetric tensor of constraints respectively. 
Let T a positive real number, Qt = Ox]0, T[ and St = <9f2x]0, T[. Consider the linear problem, 



au' + (1 - oo)Au 
a' + (w.V)<7 + udivw 
r + We{r' + (w.V)r + g(Vw, r)} 

u(0,x) 
a(0,x) 
r(0,s) 



u 



2wD[w], in Qt, 

UoO), 

<T (x), 

r (x), in 0, 
0, on St, 



(4.1) 



with 



J = F(w, 7r, ?/0 — a(w.V)w — V7r + div^, 
Q = — e _2 divw. 



and 



m 



Y < « + £ 2 7T < 29^ , in Q T . 



(4.2) 
(4.3) 

(4.4) 



4.1 Linear problem concerning the velocity u 

Consider the linear problem concerning the velocity u, 

' au' + (1 - u)Au = £, in Qt, 
< u(0, x) = u (x), in O, 

u = 0, on St, 

where Au = — (Au + Vdivu), 5 and u are given and < T < +oo. 



(4.5) 



The first Lemma concerns the existence of a unique solution of (|4,5p . By classical result of 
Agmon-Douglis-Nirenberg [I], A = —A — Vdiv is a strongly elliptic operator, and generates 
an analytic semigroup in L 2 (f2) with domain D(A) = H 2 (il) n Hq(J1) (we can see for instance 

H). 

Lemma 4.1. Let !!cl 3 o/ class C 2 , £ G L 2 (0, T; L 2 (S1)) and u € HJ(fi). Then there exists 
a unique solution of problem (f^.5| ) 

u € L 2 (0,T;H 2 (fl))nC([0,T];Hj(ft)), 
u' G L 2 (0,T;L 2 (ft)). 

Moreover, this solution satisfies the estimate 



— \\ 'II 2 

2 ll U II L 2 (0,T,L 2 (Q)) + 



^ U llL 2 (0,T,L 2 (n)) + C 1 ~~ w )H- Du l| 2 L°°(0,T,L 2 (Q)) 



+(1 -a;)||divu|| 2 LOO(0T L 2 (Q)) < 4(1 - u) \\Du \\ 2 + I 1 31 Il 2 (o,t,l 2 



(0,T,L 2 (Q)) 



(4.6) 



Proof. 

By classical result of Agmon-Douglis-Nirenberg [1], A = —A — Vdiv is a strongly elliptic 
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operator, and generates an analytic semigroup in L 2 (f2) with domain D(A) = H 2 (f2) Pi HJ(f2) 
(we can see for instance [8]). 

We start by showing the estimate (j4.6|) . Multiply ([4.5^ in L 2 (S7) by u' + a(l — u)Au, then 

f |u'| 2 + 2(l-w) f u'-Au + (l-w) 2 f \Au\ 2 
Jn Jn Jn 

= I ff-u' + Cl-w) / ff-Au. 

Integrate by parts the second term, we obtain 



|u / || 2 + (l-^)^|| J Du|| 2 +||divu|| 2 J +(l-^) 2 ||Au|| 2 < ||u'| 



+ ^||5||.||Au||. 



On the other hand, using Young's inequality on the two terms right, we get 

iisiHKll < Iw+^iiu'ir, 

(l- W )||ff||.||Au|| < i||5|| 2 + ^-^||Au|| 2 . 



Integrate over [0, T] and use the inequality 

Hdivuoll <3||L»u | 

then we get (14, 6ft . 



The second Lemma give some stronger estimates. 

Lemma 4.2 (HIE]). Under the conditions of Lemma\4l\and if dQ G C 3 , 5 6 L 2 (0, T; H 1 ^)), 
G L 2 (0, T; H _1 (fi)) and u G H^OjnHJ^). Then the solution u of problem [Jlfy given by 
Lemma \4-1\ is such that 

u G L 2 (0,T;H 3 (ft)) nC([0,T];H 2 (O) nHj(fi)), 
u' G L 2 (0,T;Hj(fi))nC([0,T];L 2 (fi)). 

and there exists a constant C 1; depend only in £1, T, a and to, such that one has the estimate 

lilt 2 , I I I |2 l I I 'I I 2 _L I I 'I I 2 

H U llL 2 (0,T,H3(n)) -T N u llL°°(0,T,H 2 (n)nHi(n)) + ll u llL 2 (0,T,Hi(n)) + ll U llL°°(0,T,L 2 (n)) 

< Ci{||AUo|[ 2 + ||5(0)|| 2 + ||3 r [|L2(o ) T,H 1 (n)) + l!^'llL 2 (0,T,H- 1 (n))}- 

(4.7) 

Proof. 

Derive in terms of t the equation f |4 . 5 j) i . then we obtain 



u" + a(l-u)Au' = &, in Q T , 
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and ujgQ(i) = for all t G [0, T]. Let v = u', then v verify the system 



v' + a(l - w)iv = in Q T , 

v(0) = v = ff(0) - a(l - a/L4u , in O, (4.8) 
v = 0, on St- 



Multiply by v the equation f|4 . 1 [) i and integrate on 0. It comes 

/ v' • v + ail —oj) Av ■ v =< v > H -i H i . 
Jn Jn ' 

After calculation, we obtain 

Id a 3a(l - w),. ||2 2 1 ,, /||2 

2 dt 4 a(l — a;) " " n 

Integrate on [0, T] and replace v and vo by their values 

III 'II 2 3a(l -u) I, , /ip 

2 I l U I lL°=(0,T,L 2 (n)) + 4 I I lL 2 (0,T,L 2 (n)) + "v 1 ^ I l Q1V U I lL 2 (0,T,L 2 (n)) 

- ^hu) W^'Ww^m + K 115(0)112 + a(1 - u) l]Au ° l 



Finally, inequality (j4.7p follows from inequality (j4.6p and (|4,9p , 

4.2 Resolution of the Transport Problems 

We consider the following two linear transport problems, 

a 1 + (w • V)<r + crdivw = — e~ 2 adivw, in Qt, 
cr(0, •) = cr ) hi 



(4.9) 



(4.10) 



and 

Jr + We(r' + (wV)r + g(Vw,r)) = 2wD[w], in Q T , 
[ r(0, ■) = t , in O, 

where a and To are, respectively, some given function and symmetric tensor defined in Q. The 
existence of solutions to this problems follows from the classical method of characteristics, (see 
for example [3, 6, [8]). The lemmas below give some estimates of the solutions of these problems. 

Lemma 4.3 (0). Let T G C\ w G L^O, T, H 3 (0)), w.n = on £ Tj and a G H 2 (ft) ; with 
o~ dx = 0. Then there exists a unique solution a G C([0, T]; H 2 (S7)) of fl^, such that 



/ cr(-,x)dx = in [0,T], 

and satisfying the following estimate 

ll°1lL°°(o,T 1 H3(fi)) - (Ikoll + ae" 2 )exp (C n ||w|| L i (0iTjH3(n)) ), 
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for some positive constant C n depending on 0. 

If, in addition, w G C([0, T]; H 2 (0)), i/ien a' G C([0,T];H 1 (J))) satisfies 



1 1 "'! lL°°(0,T,H L (fi)) - II w IIl°°(0,T,H 2 (Q)) (IKH + OLE' 2 ) exp (^C n | |w| |Li(o,T,H3(fi)) j • 

Lemma 4.4 (j3j). Lei !1 C R 3 ie a domain of class C 3 , w G I/^TjH 3 ^) n HJ(fi)) and 
r G H 2 (f2). TTiere £aere exists a unique solution r G C([0, T]; H 2 (il)) o/ fl^,il| j, suc/i i/tai 

l T llL°°(o,T,H 2 (n)) - ( ll r ol| 2 + Q^vyg) exp (p n IMIl^o/iyh 3 ^)) 

/or some positive constant C n depending on Q,. 

If, in addition, w G C([0, T]; H 2 (S1) n HJ(fl)), then t' G C([0,T];H 1 (n)) satisfies 

II t/ |!l-(0,T,H1(Q)) - C °( ll W llL-(0,T,H2(n)) + ^T^) ( ll T °ll + ^^) eXP i Cn H W HLi(0,T,H3(n)) 

4.3 Proof of Theorem I3TT1 

We are now in a position to prove the local existence of a solution to problem (|3.ip . We apply 
the Theorem of fixed-point of Schauder. 

Take T > 0, 5S l5 5S 2 > 0, and define 
<K T = {(w,n,ip), 

w G C([0, T]; H 2 (fi) n Hj(O)) n L°°(0, T, H 3 (0)), w' G C([0,T];L 2 (O)) nL 2 (0,T,Hj(fi)) 

7T G L°°(0,T,H 2 (fi)),^ G I/^O.T.H 1 ^)), 

ij) G L oo (0,T ) H 2 (Q)),V> / G L°°(0,T,H 1 (n)), 

w(0) = u , 7r(0) = o"o, -0(0) = r in O, w = in St ; 

IMlL»(0,T,H 2 («)nHj(n)) + H w llL 2 (0,T,H3(n)) + 1 1 W 1 L^(0,T,L2(fi)) + 1 1 W ' I L 2 (0,T,Hi(Q)) - 



m 

"2 



7r ll:L° o (0,T 1 H 2 (fi)) + HV ; llL°°(0,T,H 2 (n)) - ®D 

-I- I \ib'\\ < 5 

Il^o/^h 1 ^)) IlL^ccT.H^n)) — 

- < a + e 2 vr(t, x) < 2Tl 1 , in Q T }. 



Choose 25 ! such that 

> max{C 4 ||Au || 2 ,||o-o|| 2 ,||r || 2 }, (4.12) 
then (uo,o"o,t ) G £Ht- m fact, w is a solution of problem 



f w(-)GH 1 (^), 

w' + (1 — oj)Aw = 0, p.p. in 

w(0) = u , in f2, 

w = 0, on St- 

Using estimate (|4.7p . there exists a constant C4 such that 
/11 2 , 11 /11 2 



(4.13) 



l W llL 2 (0,T,H3(n)) + ll W llL°°(0,T,H 2 (n)nHi(Q)) + I l W l lL 2 (0,T,Hi(n)) + I l w l lL°°(0,T,L 2 (n)) ^ C 4 ||^4u. | 
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Thus the choose of < 3 1 in (|4,12p is enough for prove that 91t is non empty for each T > 0. 

Define now the application mapping & in this way 

A: <ft T _> X T = C([0,T];H 1 (O))xC([0,T];H 1 (^))xC([0,T];H 1 (^)) 
(w,tt,i(;) — > (u,ct,t) 

where u, a and r are solution of (|4.5p . (|4,10p and (|4,lip . respectively, with 

e 2 7r e 2 

"S = of + (1 — oj) Aw -\ Ok — wMWtt — a(w.V)w — Vtt + dWtp, 

a + e 2 7r a + s 2 tt 

Q = — e~ 2 adivw. 
If we take 

55 1 > ma X {c 4 ||Au || 2 ,e^^||(7 || 2 + ||r || 2 + l + ^-^ , C 2 (2C 5 + 1) pu (J || 2 
+C 5 pu || 4 + 3(2(1 + \\w\\ 2 c ) I |(j | | 2 + ||r || 2 ) 

+3||f(0)|| 2 + 3||f|| 2 2(0iTiHl(f , )) + 3||f|| }, (4-14) 



and 



<B 2 > e^^C 6 (||a || 2 + ||T |L + l + -^)+-^(||r |L + -^-) \, (4.15) 



2w \ 1 /,, ,, 2u 



l2 m una ■ - ■ c We ; ■ We V"-" 112 ' C,We 



and for all T small enough such that 

T - T * = min (2^(4^(1 + 1^11^)^ + 3^1' ckl) ' (4 ' 16) 

we have j?(*R T ) C SH T - 

We now use Schauder fixed point theorem. The mapping & is defined from convex, bounded 
and no empty set £Rt into Xt- To finish, we need to show the continuity of .ft in Xt- 

Lemma 4.5. To show the continuity of R in Xt, it is enough to show the continuity of A in 

2J T = C([0, T]; L 2 (S7)) x C([0, T];L 2 (Q)) x C([0, T]; L 2 (S1)). 
Proof. Let ((w n , 7i>j, VVi)^ t> e a sequence of 9^t an d tends to (w, 7T, if)), such that: 
(u n , a n , T n ) = &(w n , 7r n , i/> n ) and (u, a, r) = &(w, 7r, ^). 



Suppose that .ft is continuous in 2Jtj then the sequence |^.ft(w n , 7r n , ip n ) n j tends to K(w , tt , ip) 
in 2) T , i.e. 

lim ||(u n ,CT n ,T n )-(u,o-,r)|L =0. (4.17) 



$Ht is a compact set in Xt (see for instance [5]). Using (|4.17p . we can extract of ^(u n , o~ n , T n 
a subsequence converges in Xt to the unique accumulation point (u, a, r). Then the sequence 
((u„, a n , r n )^ = Lft(w n , 7r n , ipn)^ converges to (u, a, r) = &(w, 7T, tp) in Xt- This proved the 
continuity of .ft in Xt- D 
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Lemma 4.6. 8, is continuous in 2)t- 



Proof. 

Let ( (w n , 7r n , ip n ) ) be a sequence of an d tends to (w, n, tp), such that: 



(u n , a n , T n ) = &(w n , 7r n , ^ n ) and (u, a, r) = £(w, 7r, ip). 
Consider two systems. The first is : 



- u)Au n 


— ■Sm 




a n div w n 


= Qn, 




^w n ,r n )} 


= 2wD[w n ], 


in Q T , 


u„(0,x) 


= u (x), 




a n (Q,x) 


= o- (x), 




T n (Q,x) 


= T {x), 


in f2, 




= o, 


on St 



with 



Qn 



F(w n ,7r n , ip n ) - a(w n .V)w n - V7r n + div tp n , 
— e~ 2 div w n . 



And, the second is: 



au' + (1 - lo)Au 
a' + (w.V)er + udivw 
r + We{r' + (w.V)r + g(Vw, r)} 

u(0,x) 
a(0,x) 
r(0,x) 
u 



5, 
0, 

2wD[w], in Qt, 

u (x), 

cr (x), 

r (x), in Q, 
0, on St, 



with 



(4.18) 



(4.19) 



5 = F(w, 7r, VO — a(w.V)w — V7r + div ^, 
<5 = — e~ 2 divw. 

Let v n = u n — u, q n = a n — a and S ra = r ra — r. Using (|4.18p and (|4.19p . we obtain, in Qt: 



av' n + (1 - io)Av n 
Qn + (w„.V)g n + g n divw n 
S n + We{S' n + (w n .V)S n + g(Vw n , S n )} 



3l, 



Qi - ((w n - w).V)a - o-div (w n - w), 
Ux -We{((w n -w).V)r + g(V(w n -w),r 

(4.20) 



with the boundary conditions: 



{ v n (0,x) 
S n (0,x) 

V n 



0. 

o, 

0, in f2, 

0, on St, 



(4.21) 
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such that: 

= P(w n ,7T n ,'0 n ) -F(\V,7T,^) - a (w„.V)w n - (w.V)w - V(7T n - 7r) + div (lp n - 

Gi = -e~ 2 div (w n - w), 
%l = 2wD[w n - w]. 

First, multiply the equation (14.20^ by v n , and integrate over f2. We get: 

"(It " V "" 2 + ^ ~ W )( H Vv "l| 2 + ll divv nH 2 ) - Ikn - 7T~| Ix + IfV'n - VII? + 4 ||v n || 2 

+q 2 ||(w„.V)w„ - (w.V)w|| 2 

+ | [F( Wn , 7T n , Vn) " F(W, 7T, ^) 1 1 2 (4.22) 

We now estimate the term |[F(w n , 7r„, tp n ) — F(w, 7r, V)| | 2 on the right hand of (|4,22p . Using 
the two inequalities : 

— < a + e 2 7r{t, x) < 2Tl 1 and — < a + e 2 -ir n (t, x) < 23Jl 1} 

we obtain: 

| |F(w n , vr n , V>„) - F(w, vr, V) 1 1 2 < CVtUie 4 (1 - w) 2 I \TT n Aw n - ttAw\ I 2 + I |tt„ Vvr n - 7rVvr| | 2 

+ | [u>(7r n )V7r n — w;(7r)V7r| | 2 + | |-7r n div S n — 7rdiv S[ | 

Then, (ET221) satisfies: 



a^||v n || 2 + (l-o;)(||Vv n || 2 + ||divv n || 2 ) < C s £ n + +4 ||v n || 2 , (4.23) 



with 



|7T n -TT\\ 1 + 



- V|| 2 + a 1 1 (w„.V)w n - (w.V)w| | 2 + (1 - w) 2 | |vr n Aw n - vrAw| | 2 
+ ||7r n V7r n — 7rV-7r| | 2 + I \ w{-K n )V-K n — ii;(7r)V7r| | 2 + ||7r n div S n — 7rdiv S|| 2 . 

Second, multiply the equation (|4,20p o by e 2 q n and integrate over Q. This yields: 



e^lknll 2 < (l + C 9 |H| 2 )||w n -w|| 2 + (l + C 10 ||w n || 3 )||g r , 



< (1 + C g \\a\\ 2 ) ||w n - w|| 2 + j n \\q n \\ 2 , 



(4.24) 



with j n = 1 + C 10 ||w„|| 3 . 

Finally, multiply the equation (I4.20p ^ by S n /2oj, we obtain 



We d 
2w dt 



< 1 + 



Wed 
2uj 



2uj 11 1 



with ^ = 1 + * + ^ 



ll2 / 1 WeC 12 
|w n -w|| 1+ ^l + - + ^— ||w, 

|w n - w|| 2 + £;„ ||S n || 2 , 



(4.25) 
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The functions j n and k n , are positive and, because of the class of solutions we consider, 
j n and k n belong to L 1 (0, T). Therefore, by using (|4,23p . f[4.24|) . and (|4,25p . we deduce from 
Gronwall's lemma that: 



|2 ^ C a /■* /-4s 



K|| 2 < — / exp )e n (s)ds, (4.26) 

a Jo \ a 1 



\\q n \\ 2 < ^{dirQo + C^) exp (J j n (r)dr\ || w„(«)-w(«) | fids, (4.27) 
||S„|| 2 < (^L + CuX^ £exp(J\(r)dA ||w n (s) - w(s)|[ 2 ds. (4.28) 
The S6C[uence I (w n , 

TTri) ^n)J of IHt tends to (w, it, ip), and using (14. 26ft . t)4.27j> and t)4.28j> . we 
obtain v n , q n and S n tend to zero in 2)x- This meaning that the sequence ^(u n ,<7 n ,r„ 
m V'n) ) tends to (u, a, r) — ^(w, tt, anci is continuous in 2)x- 



4.4 Proof of Theorem ET21 

We take, as usual, the difference of two solutions (111,0*1,71) and (u 2 ,0 2 ,t 2 ) belonging to the 
class specified in the theorem 13.21 The vector function u = Uj — u 2 , the scalar function 
a = a 1 — a 2 and the tensor function r = r x — r 2 satisfy the following system: 



a 



u' + (u.V)ui + (u 2 .V)u +(l-w)iu 

<r' + (Ui.V)a + (u.V)a 2 + adiv u ± + <r 2 div u 
k r + We{r' + ( Ul .V)r + (u.V)r 2 + g(Vu l5 r) + g(Vu, r 2 )} 



with the boundary conditions: 



u(0,x) 
er(0,a;) 
r(0,x) 



u 



0, 
0, 
0. 
0. 



in 0, 
on St, 



such that: 



£ 2 - Vcr + div r, 
— e _2 div u, 
2wD[ii], 



(4.29) 



(4.30) 



d2 



F(u 1 ,cj 1 ,r 1 ) - F(u 2 ,cr 2 ,r 2 ) 



Multiply (|4.29P i . (|4.29|) o and (|4. 29j) c by u, e 2 a/a, and r/(2w), respectively, and integrate over 
f2. Summing the three obtained equations, one obtains 



1 d 



We, 



2dt 

<aC 12 



a 



2uj 



-— a||u|f + — |Mf + — |l r lf + (1 - u) ||Vu|| 2 + ||divu|| 2 +— |I T 



u, 



M + u 2 



u 



£ 

+ -c 12 



2uj 



ItTllllUilLllVull + II^HJlVullllull (4.31) 



+ \ \ a \ \ ll "!^ 1 1 cii v u 1 1 + ||(x 2 |j 2 |[(t|| |[divu|| + ||ti|| \ \a\ \ [|Vu|| + \ \a 2 \ \ \ \t\ \ ||Vu| 

We 



+ L C 12 
a 



u 



Hl3 



|cr|| + ||Vu|| ||u 2 | 



u 



1 I 1 3 



|r|| + ||Vu|| ||r 2 || 2 ||r| 
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For 5 > 0, (|4,3ip can be written as: 



~ (W + JhI 2 + ^llrll 2 ) + (1 - o,)(||Vu|| a + l|divu|| 2 )+^ ||r 




2 



(4.32) 



From (I4,32p . we then deduce that solutions (u, a, r) of ( 14, 29ft satisfy the following energy in- 
equality: 



Xs = C^dlUiH + ||u 2 || + ||u 1 || 3 )+^-(||u 1 ||^ + \\a.Wl + 2\\a 2 \\ 2 2 + \\nWl) ■ (4.34) 



The function Xg, defined in (|4,34p . is positive. Moreover, because of the class of solutions we 
consider, Xg belongs to L 1 (0, T). Therefore, choosing 5 > small enough, we deduce from 
Gronwall's lemma that u = 0, a = and r = in Qt, and that consequently u x = u 2 , cr 1 = cr 2 , 
7~i = r 2 in Qt and the system (|4,29p has a unique solution. 
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